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Abstract. ACDM, for the currently preferred cosmological density fig and cosmological 
constant predicts that the Universe expansion decelerates from early times to redshift 
z ~ 0.9 and accelerates at later times. On the contrary, the cosmological model based 
on conformal gravity predicts that the cosmic expansion has always been accelerating. To 
distinguish between these two very different cosmologies, we resort to gamma-ray bursts 
(GRBs), which have been suggested to probe the Universe expansion history at z > 1, where 
identified type la supernovae (SNe) are rare. We use the full Bayesian approach to infer the 
cosmological parameters and the additional parameters required to describe the GRB data 
available in the literature. For the first time, we use GRBs as cosmological probes without 
any prior information from other data. In addition, when we combine the GRB samples with 
SNe, our approach neatly avoids all the inconsistencies of most numerous previous methods 
that are plagued by the so-called circularity problem. In fact, when analyzed properly, current 
data are consistent with distance moduli of GRBs and SNe that can respectively be, in a 
variant of conformal gravity, ~ 15 and ~ 3 magnitudes fainter than in ACDM. Our results 
indicate that the currently available SN and GRB samples are accommodated equally well by 
both ACDM and conformal gravity and do not exclude a continuous accelerated expansion. 
We conclude that GRBs are currently far from being effective cosmological probes, as they 
are unable to distinguish between these two very different expansion histories. 
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1 Introduction 

The evidence of the accelerated expansion of the Universe coming from extensive surveys of 
high-redshift type la supernovae (SNe) in the late 1990s [1, 2] can be accurately described 
by a non-zero cosmological constant A in Friedmann equations. For reasons that might 
not be fully convincing [3-5], this simple solution does not satisfy a large fraction of the 
scientific community. Thus, more sophisticated models have been proposed: from arbitrary 
modifications of the Einstein-Hilbert action in f{R) models to brane- world cosmologies with 
extra spatial dimensions in addition to the standard four dimensional space-time, to mention 
a few [see, e.g., 6-8, for reviews]. Despite their wildly different starting point, all these 
models were conceived to describe the SN data and reproduce the expansion history of a 
standard A Cold Dark Matter (ACDM) universe. Therefore, these models are expected to 
have a decelerated period followed by the present accelerated phase. In the ACDM model, 
the transition occurs at redsfhit 1 + z = {2Q\/Qq)^^^ , where is the present average mass 
density and JIa the energy density associated to the cosmological constant. According to 
recent measures of these parameters [e.g., 9], this transition redshift is z ~ 0.9. 

Before high-z SN surveys were performed, the number of cosmological models alterna- 
tive to the standard Friedman model was very limited. Notably, the steady-state cosmology 
required an accelerated expansion more than forty years earlier [10]. Before the detection of 
the accelerated expansion, Mannheim [11] also proposed conformal gravity as a cosmological 
model, without a cosmological constant, that was alternative to the standard de Sitter solu- 
tion. After the high-z SN observations, Mannheim [12, 13] showed how conformal gravity can 
easily describe these data. Therefore, unlike recent alternative cosmologies, the steady-state 
model and conformal gravity have the remarkable property that the accelerated expansion is 
a natural feature of the model and not a requirement for building the theory. 



- 1 - 



Conformal gravity requires that the expansion of the Universe has always been accel- 
erating and a deceleration phase never occurred. Thanks to this prediction, substantially 
different from ACDM, conformal gravity is an ideal candidate to test whether the claimed 
transition between the decelerated and accelerated expansion expected in ACDM is robust. 
A potentially clean and powerful test to verify the existence of a decelerating phase in the 
early Universe is to extend the Hubble diagram of SNe to very high redshifts. Unfortunately, 
SNe are rarely observed at z larger than 1, and gamma-ray bursts (GRBs) were instead 
proposed as possible candidates to probe this high-redshift regime. The first sufficiently ac- 
curate estimates of the celestial coordinates of GRBs provided by BeppoSAX [14] enabled 
the measure of the host galaxy redshift and proved the extragalactic origin of GRBs. This 
breakthrough immediately prompted the discussion of how we can use GRBs to constrain 
the cosmological model: for example, Cohen and Piran [15] discussed how GRBs with known 
redshift could be used to estimate Oq. Atteia [16] also suggested to use GRBs as standard 
candles based on the intrinsic dispersion of various measures of the GRB brightness. 

A crucial step forward was the discovery that the luminosity of GRBs appears to be 
correlated with some of their temporal and spectral properties that can be directly observed 
on Earth [e.g., 17-19]. These correlations are not yet fully understood from first physical 
principles. However, their existence has naturally suggested the use of GRBs as distance 
indicators [20, 21]. GRBs have been observed out to redshift z = 8.2 [22, 23] and are 
expected to be observed to even larger redshifts in the future [24, 25]. They could thus be a 
powerful cosmological probe of the early history of the Universe expansion (see, e.g., [26] for 
an early review). 

The obvious procedure to use GRBs as distance indicators has two steps: (1) calibrate 
the correlations, and (2) use the calibrated correlations to estimate the GRB luminosity 
distances of a given sample. We can thus build the GRB Hubble diagram and eventually 
constrain the cosmological model. This procedure has an immediate drawback. Unlike SNe, 
there are no observed nearby GRBs and the calibration of a given correlation with a GRB 
sample requires the assumption of a cosmological model to estimate the GRB luminosities 
from their fluxes and redshifts. In principle, within the framework of a given cosmological 
model, these correlations might be used to infer the luminosities, and hence the luminosity 
distances, of other GRBs for which the remaining variable of the correlation is known. How- 
ever, using these luminosity distances to constrain a cosmological model different from the 
model used to calibrate the correlations clearly poses a problem of consistency. 

A more appropriate approach is to extract, at the same time, the correlation coefficients 
and the cosmological parameters of the model from the observed quantities. To accomplish 
this task, it is sufficient to lay out the problem within a full Bayesian context and use Markov 
Chain Monte Carlo (MCMC) simulations to compute, simultaneously, the full probability 
density functions (PDFs) of all the parameters of interest. This approach does not require 
any prior information on the cosmological model and yields results that are not plagued by 
any of the various limitations found in the literature [e.g., 27-34]. In addition, we are able to 
test cosmological models both with GRBs alone, without any priors from other probes, and 
by combining GRBs with other observables. 

Here we apply the Bayesian analysis to the GRB sample, the SN sample, and the GRB 
and SN samples combined (sect. 3), after reviewing the basic ingredients of the conformal 
gravity cosmological model in sect. 2. We discuss our results in sect. 4 and conclude in sect. 
5. 
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2 Cosmology in Conformal Gravity 



We investigate two variants of conformal gravity: the model proposed by Mannheim [11], 
that we will indicate with CG, and the kinematic conformal gravity (KCG) proposed by 
Varieschi [35]. 

CG, which is a simplified version of the original Weyl's theory [36-38], describes the 
rotation curves of disk galaxies without resorting to dark matter [39-41]. In addition, CG 
is appealing because it appears to be a renormalizable theory of gravitation [42, 43] and it 
is therefore suggestive of a possible route towards the unification of the fundamental forces; 
moreover, unlike other fourth-order derivative theories, CG does not suffer from the presence 
of ghosts [44]. However, as we mentioned above, in CG the universe has been accelerating 
at all times and never went through a deceleration phase. This very fact implies that the 
primordial nucleosynthesis lasts for a more extended time interval than in the standard model 
and the expected abundance of primordial deuterium is orders of magnitudes smaller than 
observed [45, 46]. This shortcoming of CG can be bypassed if astrophysical processes, rather 
than cosmological nucleosynthesis, supply the abundance of deuterium currently observed: 
the question is open, because the investigation of the efficiency of all the possible mechanisms 
producing deuterium is still incomplete [47] . 

KCG derives from the kinematical application of the conformal symmetry to the Uni- 
verse and implies a physical interpretation of the cosmological observables radically different 
from the standard model, as we will see below. In principle, similarly to CG, KCG does not 
require the presence of dark matter, dark energy, and an inflationary phase. In addition, 
KCG naturally explains the anomalous Pioneer acceleration [48]. 

We describe the basic properties of the CG and KCG cosmological models below. 



where C^y^x is the Weyl tensor, a is a coupling constant, and g is the determinant of the 
metric tensor g^j^^. 

The field equations are 



2.1 Conformal Gravity (CG) 

In CG^ the field action is [11] 




(2.1) 



(2.2) 



with T^'^ the energy-momentum tensor and 



6 ' 3 ' 
-R\^''' - R^".^'" + -RR'"'+ 

;A ;A 3 

-2R^^R'i + ^g'^'^Rx^R^'' - ^ 




(2.3) 



where R^'^ and R are the Ricci tensor and scalar, respectively. 
The matter action can take the form [49, 50] 




We use natural units with h — c — 1. 
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+ XS^ + i^l^r[^t^ + r^]i;-hSiJ^} (2.4) 

where S{x) is a scalar field introduced to spontaneously break the conformal symmetry, ip{x) 
is a fermion field representing all matter, r^(x) is the fermion spin connection, ^^{x) are 
the general relativistic Dirac matrices, and h and A are dimensionless coupling constants. 
The above action can be extended to include more than one scalar field [12]. The first three 
terms in the action are generated with an effective Ginzburg-Landau theory where S{x) is 
a phase transition condensate order parameter. Specifically, the term AS^ represents the 
negative minimum of the Ginzburg-Landau potential, namely the vacuum energy density 
[12]. It follows that we must have A < 0.^ 

With this action and its matter and scalar field equations, the generic energy- momentum 
tensor is 

3 6 

3 3^ 

-^s^ {r^"" - \g^''R^ - g'"'>^s^ ■ (2.5) 

Using local conformal invariance, we can set S = Sq = const and the energy-momentum 
tensor becomes 

T^" = TZ - (^R^" - Ig^^R^ - g^^XS'o (2-6) 

where = i^pj^'id" + T'']7p. 

Moreover, in the Friedmann-Robertson- Walker (FRW) metric 



ds2 = -dt' + a'{t) 



(2.7) 



1 — kr'^ 

W^'^ = and equation (2.2) implies T^'' = 0, namely 

^ (r'"' - Ig'^'R^ = TZ - ff^'^ASo^ . (2.^ 

For an energy-momentum tensor of a perfect fiuid, TZ = {p + p)u^u'^ +pg^'^, equation (2.^ 
is similar to Einstein cosmic equations recast in the form"^ 

1 1 



SvrG 



j^f^u _ ^gi^-^Rj = T^^ - gA^'^A' (2.9) 



where the usual cosmological constant A=87rGA'. It follows that, in conformal cosmology, G is 
replaced by the negative quantity —3/ (47r5Q) and the usual cosmological constant A = SvrGA' 
by the quantity — 6A5o, where A' = A5q. 

Equation (2.8) and a perfect fluid energy-momentum tensor yield the differential equa- 
tion for the scale factor a 

a2 + /t = -^-2ASoV (2.10) 



^We note that this sign is opposite to the requirement A > of Mannheim [49] and Elizondo and Yepes 
[46]. 

^We use the Weinberg [51] sign convention. 
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which is identical to the standard Friedmann equation {a? + k = SirGpa^ /3 + Aa^ /3) with the 
proper substitutions G = — 3/(47r5o) and A = — 6A5'o mentioned above. These substitutions 
make the fundamental difference with the standard Friedmann models, because in CG both G 
and A' are negative rather than positive. Moreover, G and A' depends on the same parameter 
Sq, and they depend on it in the opposite way: the lower is the gravitational cosmic repulsion 
(because G is negative), the larger is the absolute value of the vacuum energy density A'. 
Finally, CG is able to describe the flat rotation curves of spiral galaxies with the aid of a 
universal constant 70, which is related to the geometric parameter k of the FRW metric by 
the relation 7q = —Ak [40]. Therefore, we need to have fc < 0. 

In equation (2.10), we can separate the matter/energy density into a relativistic and 
a non-relativistic component: p = pnr + Pr = Pmoiao/a)^ + /3ro(«o/«)^; the dependence on 
a derives from the conservation of the energy-momentum tensor of a perfect fluid with an 
equation of state p = wp, with w = and w = 1/3 for the non-relativistic and relativistic 
components, respectively; in the FRW metric, this equation of state implies pa^^'"^^) = const. 

We define 



e. . , (2.11) 

where H = a/a. All these parameters are positive because A < and k <{). Equation (2.10) 
yields 

Oa + - enr - = 1 . (2.12) 

Unlike the standard cosmology, the expansion of the scale factor in the conformal universe 
accelerates at all times. In fact, by taking the derivative of equation (2.10), we get 

a = H^a{^ + Q, + Q!^ (2.13) 

which is always positive, whereas the deceleration parameter 

q^-'^ = -^-Qr-QK (2.14) 
2 

is always negative. 

By labelling the parameters at the present time to with the "0" subscript, we can 
rewrite equation (2.10) as 

d?a^ = H^iOAoa* + Qkoo? - O^^oa - Qro) ■ (2.15) 

Again, equation (2.15) is analogous to the standard Friedmann equation 

a^a^ = ffo(^AOQ^^ + ^feO'^^ + ^mOO' + ^ro) (2-16) 

with the standard r^nro.ro = 8-KGpnrO,ro/3HQ, Q\q = A/3Hq, and = ^nro + ^ro + ^Ao — 1- 
The main difference between the standard equation (2.16) and its conformal counterpart 
(equation 2.15) is the negative sign in front of the matter parameters 6nrO and @rO- The 
left-hand side is always positive; therefore, to have real solutions for a{t), the sum of the first 
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Figure 1. Dependence of the maximum observable redshift z^ax on the parameters A and B in 
equation (2.17). From bottom to top, the eurves are for increasing B = — A:/(pnrO*5'oao) ^^^^ range 
1 — 1000, logarithmically spaced. 

two terms on the right-hand side of equation (2.15) must always be larger than the sum of 
the latter two. On the contrary, in equation (2.16), for sufficiently small a, the matter terms 
dominate and we have the usual solutions a oc t^/^, when Om-oa is negligible in the radiation 
dominated epoch, and a oc t"^^^ in the matter dominated epoch, at somewhat larger a's. 

Therefore, in CG, a never reaches the singularity a = 0, but rather a lower limit amin > 
that is the root of the equation Oao^^ + &ko'^'^ ~ ©nro^^ — ©ro = 0. In the real Universe we 
can observe objects at very high redshift, z > 8 [e.g., 52]. Thus, we must have ciYain — 
1/(1 + Zmax) < 0.1, or possibly smaller. To obtain sufficiently small roots amim we can set 
larger and larger 6ao and smaller and smaller Bm-o and ©ro- 

To understand how amin depends on the Q conformal cosmological parameters, we can 
consider that prOi the density of the Cosmic Microwave Background, is p^o ~ 10"'^^ g cm~'^ 
[9]. To estimate pnro we consider only the contribution of the luminous component of the 
galaxies and assume a mass-to-light ratio Ad/L ~ 1, because CG, in principle, does not 
require the existence of any dark matter [39]. We find pnro ~ 10"'^^ g cm~^. Therefore 
0ro/0nrO = Pro/PrirO ~ 10~^ and we Can neglect the last term in equation (2.15), which 
becomes 

a'^a = H^QnroiAa^ + Ba-l) (2.17) 

where A = ©Ao/©nrO = —^'/PmO, and B = 0fco/0nrO = — ^/(PnrO'S'o^o)- With this approxi- 
mation, Omin = [(l + (3)^/^ + (l-Q)^/^]/(2A)^/3^ where Q = ^Jl + AB^/21A. When B ^ oo, 
f^min — ^/B, whereas when B ^ 0, cimin — 

1/^1/3. Fi gure 1 shows the dependence of z^aax 
on A and B. Clearly, sufficiently large A and B, with no specific fine tuning, can easily 
accommodate the observations of high-2 objects. 

For a sufficiently large scale factor a, we can also drop the last term in equation (2.17) 
and replace equation (2.15) with 

h'' = Hi{Qf,oa^ + eko) ■ (2.18) 
The integration of this equation trivially is 



a(t) = A/|^sinhy (2.19) 
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Figure 2. Evolution of the scale factor a{t) for different z^^x- from top to bottom the solid lines are 
for Zinax ~ 2.8, 13, 53. The dashed lines show the corresponding a{t) when the matter terms Gmo 
and Qj-o are neglected. These approximated solutions have the same HqIq of the exact solutions. 



with y = y/W\oHQt. Clearly, this solution is invalid when a — )• Omin- By imposing a(to) = 
oo = 1 for the scale factor at the present time, we have sinhyo = yOAo/Ofco and from 
q = -aajc?, we find G^o = -0Ao(l + Q'o)/q'o- 



We can determine the present age of the universe to from the relation //(to) = -f^o- We 



find 



Hoto 



atanh [\/—qo) ■ 



(2.20) 



Of course, Oao + Q/cO = because the mass parameters are negligible in this approximation. 
Fi gure 2 shows the solution ci(i) of equation (2.15) for three different .Zj^^x 

(solid lines). 

The dashed lines show equation (2.19) with the same i/o^o as the exact solutions; in other 
words, these approximated solutions have a different set of 0's but the same slope at to of 
the exact solutions. Of course, at increasing -Zmax the two solutions of equations (2.15) and 
(2.18) start to become indistinguishable at earlier and earlier times. 

In the accelerated limit, we can analytically compute the luminosity distance cIl = 
x{l + z), where the radial distance x is implicitly defined by the relation 



dr 



'o Vl — kr'^ 

and ti is the time of the light emission. One finds 



*o dt 



(2.21) 



1 



*o dt 



In 



tanh(yo/2) 
tanh(yi/2) 



whereas 



dr 



/o \/l — kr"^ 
with A: < 0. With some algebra, we obtain 



-k 



: sinh ^ { 



-kx) 



(1 + ^)^ 

qoHo 



1 + 



go 



1/2 



(2.22) 
(2.23) 

(2.24) 



Because equation (2.19) is an excellent approximation to the exact solution of equation 
(2.15) at sufficiently late times, we can safely apply equation (2.24) to the real Universe and 
consider the distance modulus m — M = ij.{z;qo) = 25 + 5logiQ[dL{z; qo) /Mpc], where go is 
the only free parameter. 

2.2 Kinematic Conformal Cosmology (KCG) 

The starting point of KCG [35] is the static Schwarzschild solution of CG which yields the 
metric 

dr^ 

ds^ = -B{r)c^dt'^ + — - + r'^dn (2.25) 
Blr) 



where 

B{r) = 1 - ~ ^^^^ - 3/37 + 7r - Kr^ ; (2.26) 



r 

f3 and 7 depend on the source mass and k is a constant. If we consider sufficiently large 
distances from the mass source [r ^ /3(2 — 3/3^)] and ignore the term f3j, that rotation 
velocities of spiral galaxies suggest to be negligible [39, 40], B{r) simplifies to 

B{r) = l+ jr - Kr^ . (2.27) 

By using the local conformal invariance, we can show that this metric is conformal to the 
standard FRW metric 

ds^ = -c^dt^ + a2(t) (^J^^ + r^dl7^ (2-28) 

with k = k/\k\ = 0,ibl, and k = —7^/4 — k. We omit the transformations between the 
coordinates (r, t) and the coordinates (r,t) which can be found in [35, 53]. 

The local conformal invariance introduces a dependence of the length and time units on 
the local metric. The redshift 

l + z = ^ 2.29 
a(t) 

can thus be interpreted as the ratio between the wavelength A(r,t) of the radiation emitted 
by the atomic transitions at the time and location of the source and the wavelength A(0,to) 
of the same atomic transitions measured on Earth now: 

Unlike the standard cosmology, where the measured redshift is due to the expansion of the 
scale factor a, in KCG, the redshift origins from the change of length and time units over the 
cosmological time and space. 

With this interpretation of the cosmic redshift, we can derive the dependence of the scale 
factor a on r or t without explicitly solving the field equations. In fact, 1+z = A(r,t)/A(0,to) 
is the ratio of two frequencies i^(0, tg) /i^(r, t) which reduces to the ratio of two time intervals, 
or the square-root of the ratio of the time-time components qqq of the metric at the two 
different locations. 

With the metric (2.25), 



l + z = J "^°°^°'^°^ = , ' (2.31) 
y -goo{r, t) ^1 + 7r- - ^r^ 
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Figure 3. Relation between redshift z and radial coordinate r of the FRW metric in KCG. The solid, 
long-dashed and short-dashed lines are for k = —1,0, and 1, respectively. The bold (thin) lines are 
for 8 = 0.06 (0.6). 

which yields, with the proper coordinate transformations to the metric (2.28), 

X + z='^ = ^V~^ -bv (2.32) 
a(r) 

where 

2 \ 1 A: = . 

Figure 3 shows the redshift 1 -|- z = a(0)/a(r) as a function of the radial coordinate r. The 
origin of r is our observer location, the coordinate r indicates locations of sources whose 
radiation has already reached us (if r > 0) or will reach us in the future (if r < 0). Clearly, 
the observed redshift z > appears only if k = —1 when r > 0. Therefore models with 
k = or 1 are not viable. The case k = — 1 also has a range of positive r where z < 0; 
spectra of sources at these locations would be blueshifted. However, this range of r decreases 
with decreasing b and can be easily accommodated in the local neighborhood of the Solar 
System. Varieschi [48, 53] actually suggests that this feature provides a solution to the 
Pioneer anomaly [54].^ When k = —1, equation (2.32) yields 



1 — 0^ 

and the two locations where z = are r = and v^^ = 2(5/(1 — 6^). There is also a minimum 
(negative) redshift at rmin = S/Vl — S^', is a real number only if |5| < 1. 

To derive the luminosity distance c^l in KCG, we consider the following argument. The 
new interpretation of redshift implies that time (i.e. l/i^) and length (i.e. A) scale as 

Al^ = (1 + z)A/o (2.34) 

At^ = {l + z)Ato , (2.35) 



*More recent investigations seem to indicate that the measured anomalous acceleration of the Pioneer 10 
and 11 spacecraft is not a gravitational effect but it rather is a thermal acceleration [e.g., 55] that might not 
even point to the Sun [56]. 



-9- 



where the subscript indicates units of the given quantity associated to objects which share 
the same location (in space and time) of the observer at the origin (namely us at r = 0) 
and the subscript z indicates the same quantities associated to objects at redshift z 7^ 0, 
measured by the same observer at the origin. For example, for an atomic transition that 
happens here we measure the frequency z^oi for the same atomic transition that happens at 
redshift z > we measure (here) the lower frequency Vz = 1^0/(1 + 2). It is important to 
emphasize that this frequency change is due to the different location in space and time of 
the atom and not to the cosmic expansion as in the standard model. The stretching of the 
space-time manifold also allows for a generic scaling of the mass: 

Am^ = /(I + z)Amo (2.36) 

where /(I + z) is an arbitrary function. Clearly, the energy IS.E oc A/^Ai~^Am scales with 
the same factor f{l + z): 

AEz = /(I + z)AEo . (2.37) 

According to these scaling laws, the relation between the luminosity Lz of a source at 
z ^ and the luminosity Lq of the same source located at z = is 

Lz = Lo^^^ . (2.38) 

In other words, Lz is the luminosity we measure on Earth when the source is at position 
r where the redshift is z 7^ 0. In practice, we do measure the flux F = Lz/'iird'j^, and not 
Lz- Classically, Lz is constant and the measured flux only depends on the distance (11- In 
KCG, Lz is not constant, but depends on the location of the source and scales according to 
equation (2.38). To account for this scaling, Varieschi generalizes the definition of flux to 
F ~ Lo/47r(f"^, where Lq is now constant, because it is the luminosity of the source at z = 0, 
and the dependence of Lz on z enters the power ay- Without a correction factor, the case 
ay ^ 2 does not of course yield a quantity with the dimensions of a fiux. Therefore, the 
proper generalization suggested by Varieschi is 

where d^s is the luminosity distance of the source at r = Vj-g where z = 0. Of course 
F{dL) = Lz/'^T:d\ still holds. By combining this relation with equations (2.39) and (2.38) 
we determine the scaling function 

In standard cosmology, the energy per unit time received on Earth is dimmed by a 
factor (1 + z)^ because of the redshift of the photon frequency and the time interval dilation. 
This dimming originates a factor (1 + z) in the luminosity distance di = aQx{l + z), where 
ao is the scale factor at the present time, and x is the radial coordinate of the FRW metric. 
In KCG, neither the energy nor the time are affected by the expansion of the scale factor. 
Therefore, the luminosity distance simply is 

dL = a(to)r (2.41) 
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where the radial coordinate r of the object at redshift z (equation 2.33) is evaluated by the 
observer at the origin, namely with b = (5(to) = b^. We thus have drs = a(to)25o/(l — <Jo) ^'^'^ 



1 + z 



25n 



5o(l + ^) + ^J(\^zY-{\-Si) 



ay 



(2.42) 



In the above equation we have used equation (2.33) with the positive sign, which is the 
solution for z > and r > r^s. 

In standard cosmology, the distance measure /x = m — M = — 2.51og]^o[-^('^L)/-^('^rcf)] = 
2.5(2 + ay) logxo(c?L/'^rcf ) clearly has ay = and d^^i arbitrarily chosen to be dref = 10 pc. 
In KCG, ay is unknown, and it is natural to choose dref = drs where z = 0. We thus obtain 



fi{z) 



2.5(2 + av)x 

5o(l + z) + V(l + z)2-(l-5o^)' 



X log 



10 



2^0 



(2.43) 



Obviously, this relation relies on the arbitrarily chosen equation (2.39) which leads to 
the definition of /(I + z) (equation 2.42). Therefore, although this choice appears to provide 
a good description of the SN data, it might not necessarily be the final correct choice [53]. 



3 Bayesian analysis 

We now apply the Bayesian analysis (see appendix A.l) to the GRB sample to derive the 
cosmological parameters of ACDM and of our two alternative models along with the coeffi- 
cients of four GRB correlations. We also derive the cosmological parameters by applying the 
Bayesian analysis to a SN sample only and to the combined GRB and SN samples. 

3.1 GRB sample 

Among the distance indicators that correlate with the GRB luminosity L = 47r(i|^Pboh where 
Pboi is the bolometric peak flux, we consider the following [21]: 

(1) . Tiag, the time delay between the soft and the hard light curve; 

(2) . trT) the minimum time over which the light curve rises by half the peak flux of the 

pulse; 

(3) . V, the variability of the light curve deflned in [21]; 

(4) . -Epcak, the photon energy where the spectral energy distribution z/F^ peaks. 

As we mention in the Introduction, we do not follow the inappropriate procedure of 
assuming a cosmological model, calibrating the relations, and deriving the GRB Hubble 
diagram. We instead apply the Bayesian analysis directly to the observables. Our GRB data 
set is {P^ov{Q'}j=i,4,z\S'}, where {Q'}j=i^4 = {Tl^g,Tl^T,V\ E'^^^^} and S' is the vector of 
uncertainties of all the measures, except for the redshift z*. 

Our task is the determination of the multi-dimensional PDF of the parameters 6 = 
{{a, 6, (Jint}j=i,4, p}, where {aj,bj) are the parameters of the four GRB correlations. In the 
ACDM and CG models 



logio -Pboi = aj + bj logio Qj - ^ogio[4.TTdl{z, p)] 



(3.1) 




0,5 1,0 1,5 3,0 0,5 1,0 1,5 2.0 0,20 0,30 0,35 0,40 0,45 



Figure 4. The marginalized PDFs of Q,q and of the parameters a, 6, and o-jnt of the L — E'pcak GRB 
correlation in the ACDM model. Black, grey, and light-grey shaded regions correspond to the 68.3, 
95.4 and 99.7 percent confidence levels, respectively. The crosses show the median values and their 
marginalized 1-cr uncertainty. The top panels show the PDFs when the Bayesian analysis is applied 
to the GRB sample alone, the bottom panels show the PDFs when the GRB and SN samples are 
combined. 

To mimic additional hidden parameters in the relations we assume that log^g -Pboi is a random 
variate with mean 

logio ^boi = aj + bj logio Qi - logio[4^4(2'> P)] (3-2) 

and variance af^^^^ (e.g., [57], [58]); finally, p is the vector of the cosmological parameters: in 
practice, p = Qq, with Oao = 1 - f^o for ACDM,^ and p = go for CG. For both ACDM 
and CG, the Hubble constant is an input parameter and we use Hq = 73 it 2 (statistical) it 
4(systematic) km Mpc~^ [60]. 

In KCG, the Hubble constant does not enter the estimate of d-L (equations 2.33 and 
2.41). The flux measured on Earth of the source at redshift z, whose luminosity measured 
on Earth is Lz = L^fil + z)/{\ + z), is 

where d^s = a(to)2(5o/(l — ^q)- Suppose the expected relation is between quantities associated 
to GRBs if they were at z = 0; in other words, suppose we expect relations of the form 

^For the sake of simplifying the algorithm implementation, for the ACDM model, we use the analytic 
approximation to dL in fiat universes provided by Pen [59]. The formula is 0.4 percent accurate when f^o is 
in the range [0.2, 1]. 
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Figure 5. Same as figure 4 for the CG model. 



log;^Q Lq = a + ftlogxo Qo- this case, we have 

logio PboiidL) = aj + bj logio Qoj - logio(4vr(i^)+ 

-»vlog,o(|). (3.4) 

Qoj is either a time tq = T2/(1 + z) or a frequency Eq = Ez{l + z):^ Tz is the time of the 
source at redshift z measured on Earth and is the frequency, measured on Earth, of the 
photons emitted at the source at redshift z. Finally, the cosmological parameters in KCG 
are p = [ay, (^o, a(to)]- 

We consider the 115 GRBs of Xiao and Schaefer [61]. The likelihood we assume for our 
Bayesian analysis is reported in appendix A.l. The analysis determines, at the same time, 
the cosmological parameters and the correlation coefficients. The top panels of figures 4, 5, 
and 6 show the marginalized PDFs of the correlation coefficients of the relation L — £'peak 
and the cosmological parameters in our three models; in KCG we show the parameter 6o. 
The PDFs of the parameters of the remaining correlations are similar. The parameters of 
the four correlations are listed in Table 1 and the cosmological parameters are listed in Table 
2. Figures 7, 8, and 9 show the four relations in the three models. 

A visual inspection of these figures and the results in the tables show that all the three 
models have no difficulty in accommodating the current GRB data. In addition, our results 

''For Ez, we use the frequency scaling relation rather than the energy relation Eg = E^/f{l + z) discussed 
in section 2.2, because the X-ray photons originated in the GRBs of our sample are detected via photoelectric 
or Compton effects, and in these electron-photon interactions the electron is ultimately sensitive to the photon 
frequency. 
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Figure 6. Same as figure 4 for the KCG model. 

for ACDM are comparable to the relations shown in [61]. Our analysis also shows that, in 
ACDM, the GRBs imply = 0.90j:o;°4 (Table 2), a factor three larger than other current 
measures. 

3.2 SN sample 

For the first time, in the section above, we have shown how the GRBs alone can be used to 
extract cosmological information without any prior information on the cosmological model 
and without any support from other observables. However, current measures of the ORB 
quantities are affected by large uncertainties and the cosmological constraints are conse- 
quently rather poor. To put more stringent constraints on the cosmological models, we can 
combine the GRBs with other cosmological probes. As an illustrative example we investigate 
the combination with SNe. Before doing so, we first apply our Bayesian analysis to the SN 
data alone. 

SNe la are not exactly standard candles. In fact, we can not trivially estimate their 
distance moduli as /.t = in"^^^ — M, where m^^^ is the apparent magnitude at the peak of 
the light curve in a given band and M an absolute magnitude valid for all SNe, because an 
empirical correlation between colour, magnitude at peak brightness, and shape of the light 
curve exists [62]. Instead, an appropriate procedure is to consider, for each i-th SN at redshift 
Zi, the following empirical distance modulus in the B band, which includes corrections due 
to the shape-parameter s of the light-curve and the color c = B — V{t = -Bmax) + 0.057 [see, 
e.g., 63, 64] 

= mf^^- - M + a{si - 1) - /3ci . (3.5) 
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Figure 7. The four GRB correlations wlien tire Bayesian analysis is applied to the GRB sample alone 
in the ACDM model. The crosses show the GRB measurements with the l-cr errorbars. The thick 
solid straight-lines are the correlations according to the coefficients listed in Table 1. The dashed lines 
show the itTint standard deviation of each relation. 



To estimate the three parameters M, a, and /3 along with the cosmological parameters 
p, we apply the Bayesian analysis to the data sample {rrii, Si, Ci, Zi, Si}, where Si is the vector 
of the uncertainties of the four observables of each SN. As in the previous section, additional 
hidden parameters are mimicked in the relation between the observable m"^^^ and fi, by 
assuming that 777,™^^ is a random variate with mean 

"^(P) = ^(P, Zi) + M- a{si - 1) + fSci (3.6) 

and variance af^^. 

Whereas for ACDM and CG, p = JIq p = qq, respectively, as in the GRB analysis, 
for KCG, p = (a, Sq) and the number of parameters is reduced by one. In the ACDM and 
CG models the expected distance modulus is = 25 + 51og((ii/Mpc), where is given by 
equation (2.24) for CG. In the KCG model, /x is given instead by equation (2.43). Therefore, 
in this latter model the value of Hq is irrelevant, whereas for ACDM and CG, /i also depends 
on the Hubble constant through di- We follow the usual procedure of removing the value of 
Hq and its uncertainty by replacing M with M' = M — 51og(-ffo)- 

The likelihood we assume for our Bayesian analysis is reported in appendix A.l. We 
determine, at the same time, the cosmological parameters, the absolute magnitude M of the 
SNe and the SN parameters a and /3, which, in turn, can be used to compute the distance 
moduli of the SNe. 



- 15 - 




We consider the 397 SNe of the Constitution set [65]. Figures 10-11 and Table 3 show our 
results. All the three models provide the same pair (a, /3) within the confidence levels; ACDM 
and CG also provide the same effective absolute magnitude M, which is fully consistent with 
the absolute magnitudes of nearby SNe [66]. Previous estimates of both Qq for ACDM (e.g., 
= 0.29 ± 0.01 by [67]) and qo for CG (go = -0.37 by [13]) are largely within 3-a of our 
estimates f^o = 0-312too23 ^^^^ 9o = — 0-225][^QQgg, respectively. We note that the reference 
value of ilo we mention was determined by combining SNe with other probes, whereas our 
estimate of Qq and Mannheim's and our estimates of qo rely on SNe alone. In addition, 
the analysis of [67], unlike ours, include the covariance matrix of each SN: based on the 
comparable results we obtain, we do not expect that our approach substantially suffers from 
this simplification. 

In KCG, the absolute magnitudes of nearby SNe also appear to be within 2-a of our 
estimate M = —16.0~^2'xi'i moreover 5o = 3.83 X 10-^ which was estimated by Varieschi [53] 
by keeping ay = 2 fixed, is within 2-a of our median value 6o = 7.81^2 8 ^ 10~^. However, 
the parameter M cannot be directly compared to the absolute magnitudes of nearby SNe 
which were estimated in a cosmological framework based on fundamentally different physical 
principles. In addition, the 68.3 percent confidence intervals are large (~ 15 percent for M 
and ^ 60 percent for 6q) because M and Jq are more strongly degenerate than M and 
or go in the other models (figure 10). At any rate, the most important result of our analysis 
of KCG, that we clearly see in the right panel of figure 11, is that the data are consistent 
with a SN Hubble diagram substantially different from what it is naively assumed to be a 
measured Hubble diagram: the distance moduli of SNe in this model are ~ 3 mag fainter 
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Figure 9. Same as figure 7 for tiic KCG model. 



than in ACDM. The reason for this result is simple: the SN apparent magnitudes are the 
direct observables, whereas the SN distance moduli are not. So, if KCG were the widely 
accepted standard cosmological model, and we tested ACDM by assuming that the distance 
moduli are the direct observables, ACDM would be completely unable to describe the data. 
Clearly our approach easily prevents us from drawing this incorrect conclusion. 

From figures 10-11 and Table 3, we conclude that the three models describe the SN 
data very well. Our result shows that SNe alone do not necessarily support an early phase of 
decelerated expansion, but are consistent with an always accelerating universe. Amendola et 
al. [68] show that the SN data are consistent with a universe where the transition between 
the decelerated and the accelerated phase can occur as early as z ~ 3. Our analysis suggests 
the more extreme conclusion that the current SN data are consistent with models where the 
transition phase never occurred. 



3.3 GRBs combined with SNe 

In the Bayesian approach, combining the SN and GRB samples simply translates into con- 
sidering a single sample S = {SNj,GRBi}, where SNj = {nii, Si,Ci, Zi}, with corresponding 
uncertainties Si, and GRBj = {P^^j, Q*-, zj, with Qj = {Tjl^g, t^T' -^peaki corre- 
sponding uncertainties Sj. As usual, our task is the determination of the multi-dimensional 
probability density distribution of the parameters 9 = {a,f3,M,af^,{a,b,aint}j,p}, where 
{a, 6, (Tint }j, j = 1)4, are the parameters of the four GRB correlations. 

The bottom panels of figures 4, 5, and 6 show the marginalized PDFs of the correlation 
coefficients of the relation L — Ep^ak and the cosmological parameters in our three models 
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Figure 10. The marginalized PDFs of the cosmological parameters of our three models and of the SN 
parameters in equation (3.5) when the Bayes analysis is applied to the SN sample alone: ACDM (upper 
row), CG (middle row), KCG (bottom row). Black, grey, and light-grey shaded regions correspond to 
the 68.3, 95.4, and 99.7 percent confidence levels, respectively. The crosses show the median values 
and their marginalized l-cr uncertainty. 



obtained when we combine the GRB and SN samples. Both the SN parameters (Table 3) and 
the GRB correlation parameters (Table 1) are barely affected by this combination. On the 
contrary, the cosmological parameters are closer to the values obtained with the SN sample 
alone than to the values obtained with the GRBs alone. Clearly, this is a consequence of the 
fact that the SN sample is larger and the uncertainties on the SN measurements smaller. As 
it happened with the SNe and the GRBs analysed separately, the three models describe the 
data properly. 

We also plot the Hubble diagrams of GRBs in our three models (figure 12), for the sake 
of clarity and to make our results more easily comparable to previous results in the literature. 
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Figure 11. SN Hubble diagrams derived with equation (3.5) and the median fit parameters in our 
three models. The error bars are computed by considering only the uncertainties on m, s, and c. 
The solid curves show the n — z relation computed with the corresponding median fit cosmological 
parameters. 

To do so, we follow the standard approach [21] and for each i-th GRB we compute /i^^ where 
j refers to one of the four relations: the directly observed quantity Q provides an estimate 
of logj^o L = a-\- 61og;^o which, in turns, returns di and thus /i^^ = 5 log]^o(^i/^Pc) + 25 
or = 2.5(2 + ay) logiQ((ii/(irs) for the KCG model. The final distance modulus is the 
weighted mean = Yl'j=i ^jA*FVSj=i where Wj = l/fj and cr| is the variance on /x^-"'^ 
derived with the usual error propagation law. Clearly, o"^. = (X]j=i 'Wj)~^^'^ is the uncertainty 
on /ij. Figure 12 shows the Hubble diagram when the parameters are derived from the analysis 
of the GRBs alone (upper panels) and when the SNe are included (lower panels). In ACDM 
and CG, the inclusion of the SNe has no visible effect. On the contrary, in the KCG model, the 
effect of the SNe is dramatic: they increase the distance moduli of GRBs by 50 percent. The 
solid curves are the distance moduli expected with the cosmological parameters derived from 
the Bayesian analysis. For ACDM our results are similar to the results of [21]. However, 
as mentioned above, both CG and KCG perform equally well, provided that the analysis 
of the data is done properly. Our result clearly shows that if we naively test alternative 
cosmological models by using GRB Hubble diagrams like those shown in figure 12 as they 
were direct measures, we can draw grossly incorrect conclusions. 

4 Discussion 

The Bayesian analysis described above shows that current data from SNe and GRBs can 
be easily described by both ACDM, that requires a deceleration/acceleration transition in 
the expansion of the Universe, and conformal gravity, that requires an always accelerating 
universe. 

To understand quantitatively which model the GRB and SN data might favour, we 
compute the Bayes factor B12 defined in equation (A. 12). We assume flat priors for all the 
parameters of the models but the internal dispersions, as explained in Appendix A.l. When 
we consider the GRB sample alone, we find that ACDM is favoured over CG by a factor 
lni?i2 = 37.9, and over KCG by a factor \nB12 = 12.0 (Table 4). When we consider the 
SN sample alone, our Bayes factor estimate shows that ACDM is favoured over CG and 
KCG by a Bayes factor lni?i2 = 6.6 and lni?i2 = 7.2, respectively. Finally, for the SN and 
GRB sample combined, the Bayes factor again shows that ACDM is favoured, by a factor 
lni?i2 = 1.5 over CG and by a factor In5i2 = 24.3 over KCG. According to Jeffreys' table 
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Table 1. Median fit parameters of the GRB relations obtained from the analysis of the GRBs alone 
(first row of each model) and from the analysis of the GRBs and SNe combined (second row of each 
model). The uncertainties are the marginalized 68.3 percent confidence intervals. The second column 
lists the number of GRBs used for the estimate of each correlation. 



[69]/ in most cases these values are large and ACDM would be "decisively" favoured over 
CG and KCG. 

However, we emphasize that this conclusion is not as robust as it might appear. There 
is in fact an extended literature on how the Bayes factor is sensitive to the choice of the priors 
[e.g., 71-76]. Investigating this issue is beyond the scope of this work, but it is clear that dif- 
ferent prior choices might substantially change our Bayes factor estimates and, consequently, 
the conclusion they suggest. 

Overall, our analysis shows that GRBs are substantially unable to distinguish between 

'^Jeffreys [69] separates the values of the Bayes factor B12 into the six intervals [(< 
1),(1,3),(3, 10),(10,30),(30, 100),(> 100)], where larger B12 values favour model Mi over model 
M2 more strongly. It became customary to indicate the intervals in natural logarithms [(< 
0), (1, 1.1), (1.1, 2.3), (2.3, 3.4), (3.4, 4.6), (> 4.6)]. The choice of these intervals clearly is arbitrary and other, 
albeit comparable, choices were also adopted [e.g., 70]. 
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Table 2. Median fit cosmological parameters of tlic GRB relations obtained from the analysis of the 
GRBs alone (first row of each model) and from the analysis of the GRBs and SNe combined (second 
row of each model). The uncertainties are the marginalized 68.3 percent confidence intervals. 
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Table 3. Median fit parameters of the SN data alone (first row of each model) and when combined 
with the GRBs (second row of each model). The uncertainties are the marginalized 68.3 percent 
confidence intervals. 
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Table 4. Values of Ini3i2, where B12 is the Bayes factor in equation (A. 12): B12 > 1 favours model 
Ml over model M2. 

two cosmological models that have very distinct predictions on the early expansion history 
of the Universe, -where the GRBs are expected to provide unambiguous results. This failure 
is mostly due to the very large uncertainties on the GRB observables. In fact, our results 
obtained by combining GRBs -with SNe confirm that the cosmological constraints are actually 
driven by the SNe rather than the GRBs [77]. Furthermore, GRB afterglo-w energies corrected 
for beaming span two orders of magnitude [78, 79], and current detectors can introduce 
significant bias against hard bursts in GRB samples [80]. It thus seems difficult to consider 
GRBs as standard candles, unless we are able to correct for these systematic effects, which 
appear to be substantial. We necessarily conclude that GRBs do not represent effective 
cosmological probes. Currently, enlarging the SN sample at 2 > 1 might remain a more 
promising method to constrain the expansion history of the Universe. 

Additional tests of conformal gravity might come from the Bayesian analysis of the 
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Figure 12. GRB Hubble diagrams based on tlie GRB data alone (upper panels) and on the sample 
of GRBs and SNe combined (lower panels). The solid curves are not fits to the data points, but the 
theoretical distance modulus of each model with the parameters derived with the Bayesian analysis. 

combination of different cosmological probes deriving from both the expansion history of the 
Universe and the properties of cosmic structures. CG seems to have difficulties in describing 
the thermodynamics of clusters of galaxies [81, 82], whereas no investigation of structure for- 
mation in KCG is available yet. However, conformal gravity remains worth investigating for 
the elegant solution it suggests to the cosmological constant, zero-point energy and quantum 
gravity problems without resorting to either dark matter or dark energy [43, 83]. 

The additional merit of the models that largely depart from ACDM, like CG and KCG, is 
that they provide predictions that are clearly distinct from the standard model expectations. 
Our results teach us that we do need a proper technique to test the data against different 
models: using a method that partly relies on the models we want to test does not generally 
yield sensible results, unless the models are very similar to each other. In this case however, 
the conclusions we draw are weak. Using a really model-independent technique returns results 
that are definitely more robust. For example, we show that KCG, where the SN distance 
moduli are ~ 3 mag fainter than in ACDM, also can describe the data. We thus confirm that 
with the Bayesian technique we can robustly and self-consistently test various cosmological 
models, no matter how different they are. 

5 Conclusion 

Current data from SNe and, more recently, from GRBs are usually interpreted as supporting 
evidence for the ACDM model, where the expansion history of the Universe has an early 
deceleration phase, down to z ~ 1, followed by the acceleration phase that lasts to the 
present time. Based on this supposedly robust interpretation of the available data, most 
alternative cosmological models were conceived to reproduce these two distinct phases. On 
the contrary, conformal gravity was proposed well before the observation of high-z SNe and 
it predicts a Universe expansion that has been always accelerating. This model is thus ideal 
to test whether the data undoubtly support ACDM. 
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We perform our test with two variants of conformal gravity: CG and KCG. KCG is 
an even more drastic alternative model than CG, because a number of physical quantities, 
including the flux from distant sources and the cosmological redshift, have a radically dif- 
ferent interpretation than in the standard model. Therefore, unlike many cosmological tests 
described in the literature, ours requires a method that does not rely on any prior information 
from the same model that has to be tested. 

We lay out the problem within a full Bayesian context, and perform three different 
analyses of a sample of GRBs, a sample of SNe and the sample of the GRBs and SNe com- 
bined. With our Bayesian approach, we simultaneously estimate the PDFs of the parameters 
of the GRB correlations, the parameters of the SN distance moduli and the cosmological 
parameters. 

Contrary to the expectation, we show that the current data can be described by CG, 
KCG, and ACDM equally well. The application of our method to the SN sample alone shows 
that the data even support the SN distance moduli derived in KCG, where their physical 
interpretation is fundamentally different from ACDM, and are ~ 3 mag fainter than in the 
standard model. Similarly, the cosmological information that can be derived from the GRBs 
does also support CG and KCG when the information is extracted properly. 

We conclude that, at face value, the ACDM expansion history is not supported by the 
data as robustly as it is naively believed, but other wildly different models can still describe 
the data satisfactorily. Therefore, current data are unable to exclude that the Universe has 
been always accelerating: both variants of conformal gravity we have investigated here are, 
from the point of view of the background expansion history of the Universe, viable alternatives 
to ACDM. 

ACDM is favoured over CG and KCG only when we resort to the Bayes factor computed 
by assuming flat priors on the model parameters. It will be crucial to investigate whether 
ACDM still remains favoured for different choices of the priors or, more interestingly, when 
we include additional cosmological probes based on the formation of large-scale structure. In 
this context, ACDM requires large amounts of dark matter that is expected to be unnecessary 
in conformal gravity. 
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A Bayesian analysis 

A.l Bayesian parameter estimation 

Consider a model M described by a set of parameters 9 with probability p{6\M) of occurring. 
The probability of measuring the set of data D when the model M is described by the 
parameters 9 is the likelihood p{D\9,M). The probability of observing a set of data D is 
thus 

p{D\M) = j p{D\9,M)p{9\M)d9 ; (A.l) 

p{D\M) is called the Bayesian evidence of the model M, p{9\M) the prior. We are interested 
in estimating the PDF of the parameters given our data set D 

For this task, we need to assume a likelihood p{D\9, M). 

For the GRBs, D = {P{^^„ {Q^.},=i,4, z\ S*}, {Q;},=i,4 = {<, t^t^ V\ i?^,,J, is the 
vector of the uncertainties of {-P^oi, {<5*}j=i,4}, 9 = {{oj, ^j, crjntj }j=i,4, p}, p is the vector of 
the cosmological parameters, and we assume the likelihood 



where 
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(A.3) 



and 



■' bol 

W; = aj+bjQ}-f{z';p) j = l,...,4 (A.5) 
where /(z*;p) is the proper function of the luminosity distance (11. W- is the mean of the 



random variate logi^Pl^ ^, whose variance is o"? , ., according to the j-th correlation. 



For the SNe, D = {mj, Si,Ci,Zi, Sj}, where Sj is the vector of the observable uncertain- 
ties. If fi{zi;p) is the distance modulus, and rrij a random variate with mean 



W, = ti{zi;p) + M -a{s,-l) + fSci (A.6) 

and variance 0";^^, the complete set of parameters is 6* = {M, a, /3, o-jnt, p}- We assume the 
likelihood 



2a' 



(A.7) 



where 



Mi 

; o'si ; cTci are the uncertainties of the observables and a^. derives from the uncertainty o"^. 
on the SN redshift Zj. 

For both the SN and GRB samples, we assume independent flat priors for all the 9 
parameters except for the internal dispersions cjint, which are positive defined. In this case 
we assume 

p(aint|M) = r^x'-^ exp(-/ix) (A.9) 
r(r) 

where x = and r(r) is the usual gamma function. This PDF describes a variate with 

mean r//Lt, and variance r/fi^. We set r = ^ = 10~^ to assure an almost flat prior. 



- 24 - 



A. 2 Bayesian model selection 

The probability of the model M to be correct, given the set of data D is, according to Bayes' 
theorem, the model posterior probability 

„M|Z,).»(^. (A.10) 

When comparing two models Mi and M2, we can compute the ratio of the posterior proba- 
bilities 

Pimp) _ p 

where 

Bu = (A.12) 
p(D\M2) ^ ' 

is the Bayes factor. If p(Afi) = p(M2), B12 > 1 clearly favours model Mi; M2 is favoured 
otherwise. Estimating the Bayesian evidence p{D\M) of a model thus provides a tool to 
compare different models for a given set of data, if the models are equally probable. 

The computation of the Bayesian evidence is not a trivial task [see, e.g., 84, for a review]. 
We use here the thermodynamic integration (or parallel tempering; [e.g., 85]). The posterior 
probability for the model parameters is 

The evidence p{D\M) is a normalization constant; we can thus define the unnormalized 
probability density 

q{e\D, M) = p{D\e, M)p{e\M) . (A.14) 

Consider the variant 

qp{e\D, M) = \p{D\d, M)fp{9\M) (A.15) 

and 

X \p(D\6,M)fp(e\M) , 
P,(.|D.M)= ^' p;(»|M) 

with the parameter /3 € [0, 1]. Because J p^(0|-D, Ad)d9 = 1, the evidence is 

Pp{D\M) = [ [p{D\e,M)fp{d\M)d9 



j qp{e\D,M)de . (A.17) 



Taking the derivative dlnpp{D\M)/dp = [l/pp{D\M)]dp,3{D\M)/dp yields 

d\npp{D\M) _ f lnp{D\e, M)qfs{e\D, M)d6 
~ Jq^{9\D,M)de 

= j In p{D\9, A4)pp{e\D, M)d9 

= {In p{D\e,M))^ (A.18) 

where the mean is intended over the set of parameters 6 and with respect to the posterior 
probability for the model parameters pi3{6\D, M). The integral jQ[dlnpf^{D\M)/d(3]d(3 = 
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lnpi{D\M) -lnpo{D\M) = \np{D\M) (because / p{9\M)de = 1) yields the logarithm of the 
Bayesian evidence entering the Bayes factor in equation (A. 12) 



A. 3 Numerical details 

For our Bayesian analysis we use the code APEMoST developed by Johannes Buchner and 
Michael Gruberbauer [86] . The first version of the code was applied to the analysis of stellar 
pulsations [87]. 

We use 2 x 10^ MCMC iterations to guarantee a fairly complete sampling of the param- 
eter space and twenty chains or values of /3 G [0, 1] . The boundaries of the parameter space 
were set to [—1000, 1000] for all the a and b parameters, [0.01, 1000] for the Uint parameters, 
[0.01,1.0] for Qo, [-1,-0.001] for qo, [1,3] for ay, [3,30] x 10"^ for 6o, and [10-^ 10] for 
a(to). The initial seed of the random number generator was set with the bash command 
GSL_RANDOM_SEED=$RANDOM. 
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